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Abstract 

It is argued that strong dynamics in the quark-gluon plasma and 
bound states of quarks and gluons is mostly due to nonperturbative ef¬ 
fects described by field correlators. The emphasis in the paper is made 
on two explicit calculations of these effects from the first principles - 
one analytic using gluelump Green’s functions and another using in¬ 
dependent lattice data on correlators. The resulting hadron spectra 
are investigated in the range T < T < 2Tc. The spectra of char- 
monia, bottomonia, light ss mesons, glueballs and quark-gluon states 
calculated numerically are in general agreement with lattice MEM 
data. The possible role of these bound states in the thermodynamics 
of quark-gluon plasma is discussed. 

1 Introduction 

The importance of nonperturbative dynamics in QCD, which is illustrated 
by the phenomena of conhnement and chiral symmetry breaking, was also 
proposed some time ago for the deconhned phase n-i In the framework 
of the Vacuum Correlator Method (VCM) (sometimes also called Stochastic 
Vacuum Model) jl] it was shown (see [Sj for a review), that four dominant 
correlators of the QCD vacuum dehne with few percent 
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accuracy all dynamics of the QCD vacuum and hadrons both below and above 
Tc, and D^^\x) is responsible for confinement, with the string tension = 
1/2/ D^^\x)(Px. Above Tc as was predicted in IB - i. D^{x) vanishes, 
while Di (x) stays nonzero and may support the nonperturbative dynamics 
at T > Tc (together with magnetic corrections due to D^^\x)^D 

At approximately the same time, starting from 1992, a careful study 
of field correlators was performed by the Pisa group pm resulting in the 
explicit forms of four independent held correlators d[^\ and 
both below and above Tc. It was concluded from these results that indeed 
vanishes at T > Tc, while three other correlators stay nonzero at least 
till T = 1.26Tc. 

It was suggested in [T] that d[^\x) is responsible for possible bound 
states of quarks and gluons in the quark-gluon plasma, which was then called 
“the strong interacting quark-gluon plasma”. Recently this phenomenon 
was observed in the lattice at T > Tc in the form of bound states of light 
qq mesons jS], heavy quarkonia (cc) 0 - ca and three-quark stetes m 
Moreover, the thermodynamic quantities associated with the QQ system, 
namely free energy and internal energy U^^^\r,T) have been 

measured o [13 cn] for the QQ distance r in the interval 0 < r < 2 fm, 
showing large asymptotic value, e.g. F{oo,Tc) ~ 600 MeV for u/ = 2 [TH] . 

These latter quantities can be explained only by nonperturbative effects, 
since perturbative OGE potential, even with increased as{r), cannot produce 
similar effect. 

At the same time, it is characteristic for the static potential Vi(r, T) 
produced by the correlator Df{x), that it indeed gives rise to the constant 
term Vi(oo, T) in the QQ interaction at large distances, which can be viewed 
upon as the sum of constant selfenergies of Q and Q. 

It was argued in the recent paper HZI by one of the authors (Yu.S.), that 
this strong interaction observed in the quark-gluon plasma on the lattice in 
i-ra. and possibly seen in the ion-ion collisions at RHIC can be ex¬ 
plained by the correlator Df{x). To check this prediction in na the analytic 
form of Tf (x) was used calculated before in m, exploiting the connection 
of held correlators to the gluelump Green’s functions. In this way the static 
QQ potential Vi_(r, T) was computed from Tf (x) and compared to the lattice 
free energy F^^^\r,T). The resulting agreement has demonstrated that the 
D\ (x) can be considered as a calculable source of strong interaction above 
Tc. 
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Based on Vi(r, T) the possibility of bound states of heavy quarkonia, 
glueballs and baryons was established, using the Bargmann condition for the 
static potential. 

It is a purpose of the present paper to calculate explicitly the possible 
bound states of quarks and gluons using directly the previously calculated 
d[^\x) on the lattice as an input. In doing so a new £t of lattice data is 
done, using the analytically calculated form of d[^\x), found in [HlEj. For 
comparison we also use static potential obtained from the analytic expression 
of d[^\ found in [1^] for T <Tc and analytically continued to T > Tc. 

With the help of these sets of static potential, lattice and analytic, 
and ^ -yve calculate the bound states of charmonia, bottomonia, strange 

quarkonia and glueballs and compare it with existing lattice data. 

Having calculated possible bound states, we discuss their possible role in 
the thermodynamics, and in particular, the effective masses of quarks and 
gluons and the contribution to the free and internal energy with the aim to 
explain the difference (e — 3p) above Tc and ratio e/esb- 

The paper is organized as follows. In the next section we introduce the 
analytic gluelump expression for D\ and the static potential derived from 
it. In section 3 the analytic form of is tested with existing lattice data 
in the interval l.OOTTc < T < 1.26Tc and in section 4 the static potential is 
derived from the lattice data. In section 5 we write the Hamiltonian for heavy 
and light quarkonia, glueballs and quark-gluon states. In section 6 numerical 
results for bound states in concrete systems are considered and discussed in 
comparison to lattice data. The possible role of calculated objects in the 
thermodynamics of quark-gluon plasma is envisaged, and the suppression 
factor for the colored states is derived. 


2 Analytic form of field correlators and the 
QQ static potential at T = 0 

We start with the standard dehnitions of the held correlator functions D{x) 
and Di{x), dehned as in |1], which are measured in 

-j^{trf{F^^{x)^{x,y)Fxa{y)^{y,x))) = D^^^xa{x,y) = 
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= - 6^„5^x)D{x -y) + ^ + perm^ Di{x - y) ( 1 ) 


where h^ = x^— y^, $(x, y) = Pexpig fj A^dz^ and trj is the trace in the 
fnndamental representation. Our hnal aim in this section will be to connect 
D{x), Di{x) to the gluelump Green’s functions. These Green’s functions are 
however not accessible for direct analytic calculation and to proceed one 
needs to use Background Field Formalism (BFF) [20], where the notions of 
valence gluon held and background held are introduced, so that total 
gluonic held is written as 


Ajj, — ( 2 ) 

Using (ED one can write the total held operator F^y{x) as follows 
F^,y{x) = d^A^ - dyA^ - ig[A^, A^] = 


T Bi,^ T ig^Qi^ T (Xy T By^ (3) 

= bf,ay - byttf, - ig[a^,, ay] + 

Here the term, Pjiy^ contains only the held To proceed one can use the 

idea j2I] of the background held B^!’^ as the collective held of all gluons with 
color indices b, with b occupying most of indices from 1 to — 1, while the 
color index ”a” of a“ having only few values. The physical idea, realized in 

naini is that after averaging over all helds Bj!’^ in the vacuum, the resulting 
interaction for the valence gluon is diagonal in index ” a” and has the form of 
the white adjoint string. It is clear that when one averages over held a“ and 
sums hnally over all color indices a, one actually exploits all the helds with 
color indices from so that the term can be omitted, if summing 

over all a is presumed to be done at the end of calculation. In this section 
we shall concentrate on the hrst two terms on the r.h.s. of ©, which yield 
Di{x). 

Assuming the background Feynman gauge, = 0, we shall dehne now 
the gluelump Green’s function as 

G^,y{x,y) = {traa^,{x)^{x,y)ay{y)). (4) 
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We can also write for the gluelump Green’s fnnction ( 0 )) 
— uY) in the limit of vanishing glnon spin-dependent 
interaction, see discnssion below. 

As a resnlt one obtains from dS© the following connection of and 
/((x — nY) (see [in] for details of derivation) 


D 


( 0 ) 

42,4Ai 





d d 
dx4 dy^ 




yr) + 


d d 

dxi dyk 





(5) 


on the other hand nsing ( 1 ) with = x^ — y^ one can express throngh 

Di as 

<«('») = Wh) + i ( 6 ) 

and for hj = 0, i = 1, 2, 3, ^4 7 ^ 0 one obtains 

Di{x) = - ^2 ^ G^.u{x,y) = 6f,J{{x-yY) ( 7 ) 

To obtain information abont the glnelnmp Green’s fnnction G^i, one 
can nse the path-integral representation of G^Y^,y) in the Fock-Feynman- 
Schwinger (FFS) formalism (see |22| for reviews and original references), 
which was exploited for glnelnmp Green’s fnnction in [22j 

poo 

Gf,u{x,y) = tVa ds{Dz)^ye~^ {Wl£\G:,y)) ( 8 ) 

Jo 

where K = ^ dr and 


^yu^Gxy) = PPpl^y^^iiig j Axdzx)exp{2ig drF^p( 2 :(r))) | (9) 

and the closed contonr Gxy is formed by the straight line from y to x dne to 
the heavy adjoint sonrce Green’s fnnction and the path of the valence glnon 
from x to y. Note that the nontrivial {yv} dependence of the r.h.s. of 
o occnrs only dne to the F^p] expanding in powers of this term, one has 
Wj = + ... In what follows we shall neglect W^) with 

n > 1 since these terms correspond to the effect of spin-dependent forces in 
glnelnmp, which are relatively small and were acconnted for in |23j . 
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Neglecting in gluon fields altogether we obtain the perturbative re¬ 
sult, ^ 




At^‘^{x — yy 


( 10 ) 


This is the leading term in the expansion of at small \x — y\, while the 
next order term is found in na to be 


D^x) = 


4 ( 720(5 


TT 


1 TT^Go 


x^ 


+ 


2ANr 


+ ... 


where G 2 is the standard gluonic condensate pa 


D(o) + r>i(o) = 




( 11 ) 


( 12 ) 


One can prove consistency of the resulting Di{x) as follows. 

Checking constant term in (HH), one can compare T*i(0) on the l.h.s. 
of (HH) with the r.h.s., T)i(O) = ^ ■ ^G 2 = ^(T>(0) + T>i(0)), where 
Zi)(0) + T*i(0) on the r.h.s. of (|TT| are defined by the gluon condensate, Eq. 
m- Since < 1, this estimate of T*i(0) is reasonable and suggests that 
for as = 0.4 the magnitude of T*i(0) is 0.2 T>(0). 

This ratio is in agreement with the lattice calculations in 
In another check one considers the singular term, in_ 

serts it in the static QQ potential. The static QQ potential can be expressed 
through D and Di, as was done in [23]: 


V{r) = 2r / dX \ dvD{X,v) + / XdX 

Jq Jq Jo 

= VD{r) + Vi{r). 


di'[—2D{\, u) + Di{\, u)] 

(13) 


Inserting in ca the perturbative part of Di from cu one obtains the stan¬ 
dard color Coulomb potential Vc(r) = —thus checking the correct nor¬ 
malization of Di{x). 

Another form of Gik{x, y) is available at all distances and practically 
important at large \x — y\, namely 


Gid^.v) = NdN^ - 1) J]*«>(0)<'=r(0)e-""l“'->'l (14) 

n=0 
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where are eigenfunction and eigenvalue of the gluelump Hamil¬ 

tonian, [2H], details are given in Appendices 1,2,3,4 of ref. ng. 

For one can use the known equation, which is obtained from the 

eigenfunctions of Hq through the connection 

( 15 ) 

Inserting m into m one obtains 


n=0 


{(^adj) M„\x-y\ 

471 


(16) 


It is clear that for x ^ y the sum in P|) diverges and one should use instead 
of (jl6j) the perturbative answer (uni). For large \x — y\ one can keep in (IT?)|l 
only the terms with the lowest mass, i.e. for the color electric gluelump state 
1 , which obtains for spacial = i,k 

Thus one gets 

||.-,Hoo ^ (iVc - (17) 

The eigenvalue Mq was found in j2Sl to be Mq = (1.49 0.98) GeV for 

(7/ = 0.18 GeV^ and = (0 -h 0.195). This should be compared with the 
value Mo = 1.5 ± 0.4 GeV, obtained from the QGD sum rules in [2^] and 
with the lattice value Mq 1 GeV obtained in iniizi from asymptotic of Di. 

Using m one can dehne from (HZl) the nonperturbative part of Di, which 
is valid at large |a:|. 


D 


(nonp) 

1 


(x) 


G2{f)ois‘^Moaadj G2{f) 


2N, 


and the total Di due to dZD and (uni) can be represented as 


( 18 ) 


D^{x) = + Z)l"‘'“»>(i) (19) 

TTX 

and 7 plays the role of screening length of gluon. 

One can easily see that insertion of m into Vi{R) allows to express Vi{R) 
in terms of /(x), 

2 poo 

Vi(R) = -jpJ^ + (20) 
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For the purely perturbative /(x^), 



( 21 ) 


one has from (HI the standard lowest order Coulomb interaction 



( 22 ) 


Note that the last term in this case, Fi(cxo) = is diverging 

at small v, which is connected to the well-known perimeter divergence of 


hxed-contour Wilson or Polyakov lines, which is renormalized on the lattice 
[121OE], subtracting short-distance part of potential at T = 0. In what follows 
we shall concentrate on the nonperturbative part of Vi{R) representing it as 


Vi{R) = where I//’®’’* can be treated as the screened 


Coulomb potential, so that renormalization of the Wilson (Polyakov) lines 


would affect only 


One can easily see in (EOl) that Pi(cxo) is actually the sum of self-energy 
parts of quark and antiquark, occurring due to the gluelump exchange with 
time interval v. 

3 Analytic form vs lattice data at T > Tc 

A detailed study by numerical simulations on a lattice of the behaviour of the 
gauge-invariant two-point correlation functions of the gauge-held strengths 
across the deconhnement phase transition (T ~ T,,), both for the pure-gauge 
SU{3) theory and for full QCD with two havours, has been performed in 
Ref. j2Z] • Quenched data published in j2Zj agree within errors with previous 
determinations j2Hl (obtained on a 16^ x 4 lattice, in a range of distances 
from 0.4 to 1 fm approximately) but have been obtained on a larger lattice 
(32^ X 6 ) and with much higher statistics. 

For the beneht of the reader (and essentially for dehning the notations 
used in the rest of the paper) we report here some technical details about 
the lattice determination of the correlators in M- To simulate the system 
at hnite temperature, a lattice is used of spatial extent 3> W, being 
the temporal extent, with periodic boundary conditions for gluons and an- 
tiperiodic boundary conditions for fermions in the temporal direction. The 




temperature T corresponding to a given value of /3 = 2Nc/g'^ is given by 

TV, ■ a = ^ , (23) 

where a is the lattice spacing. In the quenched case a only depends on the 
coupling (3 and, from renormalization group arguments, 

o(/3) = . (24) 


where f{(3) is the so-called scaling function and is the scale parameter of 
QCD in the lattice regularization scheme. At large enough /3, /(/3) is given 
by the usual two-loop expression: 



51/121 

exp 


[1 + 0 ( 1 / 1 ()| 


(25) 


for gauge group SU{3) and in the absence of quarks. The expression 
can also be used in a small enough interval of /3’s lower than the asymptotic 
scaling region, and then A^, is an effective scale depending on the position 
of the interval considered. The lattice used in Ref. izzi for the quenched 
case was a 32^ x 6 (in our notation, = 32 and A, = 6) and the critical 
temperature T^. for such a lattice corresponds to f3c — 5.8938 [221 • The range 
of values of /T’s considered in Ref. Ea goes from [3 = 5.85 to {3 = 6.10 and 
in this interval the effective scale A^ is about 4.9 MeV. 

At hnite temperature (Ao- 3> A,) the 0(4) space-time symmetry is bro¬ 
ken down to the spatial 0(3) symmetry and the bilocal correlators are ex¬ 
pressed in terms of hve independent functions n-i. 00 Two of them are 
needed to describe the electric-electric correlations: 

y)Ek{y)<!>{y, a:)]) 


= h 


ik 


D^ + Df + ul 


dDf 

dul j 




dDf 

dv? 


(26) 


where Ei = E^ is the electric held operator and = x^ — y^, [if = [x — y)"^]. 
Two further functions are needed for the magnetic-magnetic correlations: 
-,2 
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— {trf[Hi{x)^{x, y)Hk{y)^{y, x)]) 


A 


= h 


ik 


'iH I I ^2^Ei 


o" + Of + r 


dv? 


XliXlk 


aof 

du^ 


(27) 
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where Hk = ^SijkFij is the magnetic held operator. 

Finally, one more function is necessary to describe the mixed electric- 
magnetic correlations: 

^{trf[Ei{x)^{x,y)Hk{y)^{y,x)]) = --SiknUn ■ (28) 


In Eqs. dSni), pTjl and (j28|l . the hve quantities ^ Df, and Df® 

are all functions of iP, due to rotational invariance, and of m|, due to time- 
reversal invariance. 

The following four quantities have been determined in 


Df(u^0) = D^{u‘^,0) + Df{u\0) + u^^{u^,0) ; 

Df(u2,0) = D^{u‘^,0) + Df{u\0) ; (29) 

Df(M^O) = T>^(m^0) + T>f(M^0) +m^^|-(m^0) ; 

Df(M^O) = T)^(M^0)+ T)f(M^0) , (30) 


by measuring appropriate linear superpositions of the correlators (i2nD and 
at equal times {u^ = 0). Concerning the mixed electric-magnetic cor¬ 
relator of Eq. (I28j) . it vanishes both at zero temperature and at hnite tem¬ 
perature, when computed at equal times = 0), as a consequence of the 
invariance of the theory under time reversal. 

The results found in Ref. EH, both for the quenched and the full-QCD 
case, are in agreement with those already found in Ref. I2HI and can be 
summarized as follows: 


(1) In the conhned phase (T < Tc), up to temperatures very near to Tc, 
the correlators, both the electric-electric type (Ei and the magnetic- 
magnetic type (EH), are nearly equal to the correlators at zero tem¬ 
perature: in other words, ~ ~ D and Df ~ ~ Di for 

T < Tc. 

^The definition of the gauge-invariant field-strength correlation function 
adopted in m differs from the one given in this paper by the absence of the multi¬ 
plicative factor 1/A^c on the left-hand side of Eq. (^. Therefore all functions , Df,... 
used in this paper are smaller then the corresponding functions used in m by a factor 
iVc = 3. 
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(2) Immediately above T^, the electric-electric correlators (I29|l have a clear 
drop, while the magnetic-magnetic correlators PI|) stay almost nn- 
changed, or show a slight increase. 

(For the quenched theory the behavionr of and Df is shown in Figs. 1 
and 2 of Ref. [2Z1 respectively, at different valnes of T/Tc with the physical 
distance in the range from ~ 0.25 fm np to ~ 1.25 fm. The analogons 
behavionr for and is shown in Figs. 3 and 4 of Ref. EH.) 

Moreover, a best-£t analysis of the data has been performed in [2Zj , both 
for the quenched and the fnll-QCD case, with fnnctions for D and Di having 
a pertnrbative term a/x‘^ (here x = Ini) plus an exponential non-perturbative 
term Aexp{—iix): this analysis has shown explicitly that the electric gluon 
condensate drops to zero at the deconhning phase transition. 

In this section, inspired by the results found in and partially reported 
in the previous sections, see Eqs. m and (O, we shall present the results 
obtained by performing alternative best hts to the quenched lattice data for 
the electric correlators (P|) at temperatures T above the critical temperature 
Tc (and at equal times, i.e., M4 = 0), with the functions: 

D^ix) = 4 ’ (^) = -e-^^ + 4 , (31) 

X^ X x^ 

where, of course, all the coefficients must be considered as functions of the 
physical temperature T. According to the results found in Refs. (271I2HI, the 
function is taken to be purely perturbative (so behaving as l/x^) in the 
deconfined phase (T > Tc). Indeed, from the conclusions of Refs. [H 12 E) , 
one expects that the non-perturbative part of is related to the (tem¬ 
poral) string tension and should have a drop just above the deconfinement 
critical temperature T^. In other words, the non-perturbative part of is 
expected to be a kind of order parameter for confinement and this is fully 
confirmed by the results found in Refs. I2III2H1 On the contrary, Df does 
not contribute to the area law of the temporal Wilson loop and we use for 
this function a parametrization derived from Eqs. (HHD and (d, consisting 
in a sum of a non-perturbative term B/x exp(—Ma:) plus a perturbative 
term h/x^. (As already noticed in |27]) the perturbative coefficients a and 
b are regularization-scheme dependent. In Eq. m we refer to the lattice 
regularization scheme; other schemes could give different values ISO]- We will 
comment again on this question in the next section.) 

The best fit has been performed to the data for both (perpendicular 
and parallel) electric correlators (I29|l for distances from 3 up to 6-7 lattice 
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spacings (corresponding approximately to the range of physical distances 
0.3 -h 0.7 fm), i.e., for those distances where we have data for both parallel 
and perpendicular electric correlators at all temperatures (see Figs. 1 and 2 
in Ref. |77ji. 

After having tried many hts with all the parameters free for each given 
temperature T, which were rather unstable and not strictly conclusive, we 
have taken the mass M of the non-perturbative term of Df in m and 
the perturbative coefficients a, b to be temperature independent (at least in 
the range of temperatures that we are considering): indeed, this fact, i.e., 
the temperature independence (in the short range of T considered) of the 
perturbative coefficients and of the correlation length of the non-perturbative 
terms, was hrst suggested and conhrmed in Ref. EH." 

In Table I we report the results obtained by htting simultaneously all 
the data for the perpendicular electric correlator Df [see Eq. d^ ] at tem¬ 
peratures T > Tc with the functions (EH), where only the non-perturbative 
coefficient B is considered to be temperature dependent. The value of the 
mass M comes out to be about 1 GeV and the non-perturbative coefficient 
B drops rapidly to zero (within the errors) going from T/T^ = 1.007 to 
T/Tc = 1.261. The value of the perturbative coefficient a -|- 6 is perfectly 
consistent (within the errors) with the value 0.90(3)/Ac found in ^7\. In 
Table II we report the results obtained by a best £t similar to the previous 
one, but fixing the value of the perturbative coefficient a -|- 6 to the value 
0.90/A7 = 0.30 found in [77] . 

Finally, we have performed a best £t to all the values for the difference 
Df{x) - Df{x) = (32) 

between the two electric correlators (|2ni) at T > Tc. The results are reported 
in Table III and, within the very large errors, they roughly agree with those 
obtained in the two previous best hts. Let us observe, in particular, that 
the value of the perturbative coefficient b agrees, within the errors, with the 
value 0.35(1)/A^c found in [21] • If we repeat the best £t by hxing the value of 
the perturbative coefficient b to the value 0.35/A"c ~ 0.12 found in |27|, we 
obtained the results reported in Table IV. 

^We remind again, however, that the non-perturbative terms for the functions D and 
Di in Ref. |27| were taken to be exponentials, Aexp(—/ri), with the same correlation 
length A = 1/^ for all functions. 
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The values of for the best fits considered are satisfactory and we can 
conclude that the functions (ED), inspired by Eqs. m and (uni), represent a 
reasonable parametrization of the correlators in the deconfined phase. 


4 Static potential from lattice correlator 
measurements 


In this section we use the results from the best fits to the lattice data for 
the electric correlators (Ei above with the functions m, that we have 
obtained in the previous section, for deriving the QQ static potential Vi{R, T) 
in the deconfined phase. 

According to the results found in na, and partially reported in the pre¬ 
vious sections, the expressions m for the electric correlation functions 
and Df imply for the perturbative and the non-perturbative part of the 
static potential eA T > Tc, 

fo {R, T) = (i?, T) + (i?, T) , (33) 


the following approximate expressions: 


V(P^rt)^R^T) 


Trb 

Ir 


1 


—arctan(i?T)-In 

vr TT 



(34) 


(satisfying the condition = 0) and: 




B 


Ki{MR) MR - + MR) 




B 


T 




(35) 


where Ki{x) is the modified Bessel function. 

As already observed in the previous section, standard perturbation theory 
provides the following estimate for the coefficient b [compare the parametriza¬ 
tion m for Df with the expressions m and CT]: 


b 


4a.C2(/) 

TT 


+ 


C2{f) 


2iV, 


(36) 
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and the first (zero-temperature) term in the right-hand side of (Id4|l . i.e., 
= 0) = —7rb/{AR), is nothing but the standard Coulomb poten¬ 
tial written in Eq. (El- 

When extracting the parameter b from lattice measurements of the cor¬ 
relators, as we have done in the previous section, the coupling constant 
in dSHD must be identihed with the bare lattice coupling constant, which, ex¬ 
pressed in terms of (3 = 2Nc/g‘^, is given by Nc/{2tt(3). In our case (W = 3) 
the range of (3 values goes from 5.9 to 6.1, corresponding to an ~ 0.08. 
Using this value of as, one immediately verihes that the values of b reported 
in Tables III and IV agree almost perfectly with the perturbative estimate 

fjHUj) . 

In Fig. 1 we show the behaviour obtained for the static potential Vi{R, T), 
given by Eqs. (Ei-dnni), as a function of the distance R, for different values 
of the temperature T > T^, using for B{T) and M the central values reported 
in Table I and for the parameter b the perturbative estimate dsni), with the 
bare lattice coupling constant as = 3/(27r/3), as previously discussed. 

The error associated with each curve in Fig. 1 can be estimated by 
observing that the large-distance behaviour of the potential is dominated 
by the non-perturbative part u/"'^^(cx)) in Eq. given approximately by 
where, on the basis of the results reported in Table I, the value of 
M is known with a relative error of 10% and the value of B (for each T) 
is known with a relative error of at least 30%. As a consequence, the non- 
perturbative part of the potential can be derived from the results of Table I 
with a relative error of at least: 


^V^np) 


~ 50% . 


(37) 


5 Equations for bound states 

In this section we derive Hamiltonian for different quark-gluon systems, using 
static potentials (d, derived for the color singlet QQ system, which for hnite 
temperature T has the form nn, 

^1/T j^r 

U/^«V,T)=/ duil-uT) (38) 

Jo Jo 

We can now generalize the static interaction to the case, when two color ob¬ 
ject A and B with the Casimir coefficients C{A), C{B) combine into common 
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color state D, with Casimir C{D). The answer is^ 




1 

2 ^ 


{C{D)v}^^\oo,T) + (C'(A) + 


+C{B) - C{D))vl^^\r,T)} = aVi^^'>oo,T) + hVl^^\r,T). (39) 

Here C{f) = | is the Casimir number for the fundamental charges. To 
illustrate (Ei several examples are considered in Table 5. 

For multicomponent systems one can similarly hnd, 

°V'M»c>(r,,r2,r3,T) = -^{C(D)vi'^^\oo,T)+ 

+ {Ca + Cb + Cc-Cd)\Y. - ry. T)} («) 

i>j 

in particular 

(QQQ)i : = i y«>(ri - r,. T) (41) 

i>j 

(QQQU -■ ho®®®’ = - 1 - rj,T) (42) 

i>j 


{QQQ)t ■ h®®®’ = sh®®’(aj.r) 


5 ;h«®'(ri-r^,r) 


(43) 


i>j 


One can check that at large distances all systems consisting of uq quarks 
or antiquarks and Ug gluons tend to the constant limit, independent of D 


D 


Vo(noO.»,»)(|r. _ r ^ + E^n„ 


SQ = ih''®’(«3,T), B, = 5h«®>(oo,r). (44) 

■^The construction is similar to that found in PSl for the coefficient b of Vi^^\r,T) 
but differs in the presence of the first term proportional to a, since in m constant term 
was not taken into account. 
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Since nonperturbative part of T) ~ 0{r‘^),r —> 0 one obtains 

the lower bound on the nonperturbative part of which is 

g£iv'/o®(oo,r). 

As a consequence, one can predict the absence of bound states in some 
channels, e.g. in {QQ)s, {QQ)^, {QQQ)io etc. 

As one application of the general relation iOD we show in Fig.2 the static 
potentials r 2 , rs, T) = T) of three static fundamental quarks 

in three different representations D: singlet {C{D) = 0), octet(C'(Zi)) = 3) 
and decuplet, {C{D) = 6) in the symmetric conhguration with Vi = R,i = 
1, 2, 3, as function of i?.One can see in Fig.2 that all three potentials tend to 


the same limit at large R, in accordance with Eq. (jUj), while 

deviations from asymptotic at all distances are proportional to for 

singlet, decuplet and octet respectively, as prescribed by Eqs. fl4HI43|) . This 
is in agreement with lattice calculations of free energies Fgqq\R, T) presented 

in (32]. 


Having constructed static potentials for different systems, we can now 
exploit the relativistic Hamiltonian technic, developed in j2S| and successfully 
used for mesons, baryons, glueballs and hybrids in the conhnement phase (see 
[23j for a review). This technic does not take into account chiral degrees of 
freedom and is applicable when spin-dependent interaction can be treated as 
perturbation. Therefore below we stick to the Hamiltonian technic of IHHI and 
consider heavy quarkonia, and baryons, leaving light quarkonia with chiral 
symmetry restoration to another publication. 

Leaving details of derivation to |23IS11, one can write the bound-state 
equation as 


H'ipri ^n^ni 


(nq+rig x 2 .\ _ ) 

= mm I ^ ^ ^ -RnQ+ng) + aV^^^\(X),T) | , (45) 

where we have introduced the einbein variables /Xj for quarks and gluons 
with the stationary values playing the role of the constituent masses. 
The Hamiltonian has the form 

H = Ho+ Hs + Hse (46) 
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where 


Ho= + ( 47 ) 

i=i 

and Hs, Hse are spin-dependent and self-energy parts of Hamiltonian defined 
in terms of field correlators. 

The search for bonnd states in the QQ and QQQ systems was done using 
two types of interaction v}^^\r,T). 

i) The hrst one based on the analytic representation (UHl), (HI obtained in 
the conhned region, and analytically continued into deconhnement as 

VQQ(r, T) = T) + Vf{r, T) (48) 

with 

pl/T 

v}'^P\r,T) = a{T) du{l - uT)[e''^^ - (49) 

Jo 


Aa 


l/i^(r,T) = -^(p(r,T), «, = 0.3 
with the temperature-modihed Coulomb term as in (ED 

1 


2 tT 

ip{r, T) = ( 1-arctan(rT)-In 

TT TT 


1 + 


(rT)2j 


„-7r 


(60) 


(51) 


Here the coefficient a(T) is 

a{t) = ao- 0.36^—^, (52) 

^ C 

and the value of Oq coincides with that obtained in dUD for the con¬ 
finement region, 


ao = 2C2{f)asaadj = 0.648 GeV^ (53) 

while Mq is taken in one case the same as the lowest gluelump mass 
Mq = 1 GeV |2S1, and in another it was assumed to be decreased in 
the deconfined region to the value Mq = 0.69 GeV. The value of 7 was 
taken at 7 = 0.2 GeV and 0.69 GeV. 
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ii) The second choice is based on the lattice determination of the correla¬ 
tor Df{x), done in 1231211 and described above in sections 3 and 4, 
and in Fig. 1. Here one must rescale the resulting Df{x), taking into 
account that lattice coupling a* ~ 0.08. From theoretical dehnition 
of through the gluelump Green’s function in m one can 

see, that is proportional to (with higher order terms propor¬ 

tional to a”, n >2 times more complex gluelump Green’s functions, as 
discussed in m)- Since we are interested in large distance behaviour 
of T), one should use the infrared-saturated value of ^^(cxd), 

which in the confined region (and hopefully sX T > Tc but close to 
Tc) can be taken from analysis of meson spectra j2S] and background 
perturbation theory (20] • According to this the acceptable value is 
as(cx)) ~ 0.6 -h 0.45 for nj = 5 -=- 0. 


In Fig.l a rescaling parameter of the lattice defined Di 


{nonp) 


X 


and 


is taken as .^ = 


Ctsjcxj) 


as (lattice) 


0.6 

0.08 


= 7.5. 


In Fig.l on the r.h.s. of the hgure, the rescaled values of are given, 
and one can hnd the maximum value of v}'^^\oo,T) at T = 1.007 is 
equal to ~ 0.53 GeV. These values are in farly good agreement (within 
the large errors) with the lattice measurements of Fi{r,T) for T > 
reported in HE]. 


The two interactions, described in i) and ii) respectively, and denoted 
Vi{I) and Vi(//), have been used to calculate the bound states of different 
binary systems, both white and colored, as shown in Table 6. 

The mass of the binary state can be computed according to (I45M7|1 . with 
inclusion of the total selfenergy for colored bound states. 


Mi = aiV^'^^\oc,T) + 


m 


1 + + /^i +/^2 ^ ^(i) 


\2/ii ' 2/i2 ' 2 

where the eigenvalue (/i) is a solution of equation 


+ eW(^) (54) 


^ T) ) i,{r) = e»(A)V'(r). 


( 66 ) 


Here ft = and mi, m 2 , hj, bi are listed in Table 6. 

Several words about the choice of masses mi, m 2 in Table 6. In the 
conhned phase at T = 0, mi, m 2 are current masses and their values for charm 
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and bottom quark are assumed to be the same for T > and correspond 
to the values listed in PDG |2S| and used for heavy quarkonium spectra in 
jd7j . For light quarks and gluons for T > one should take the values either 
found from the lattice analysis [3H] , or from the quasiparticle calculations of 
quark-gluon plasma [ 2^1 (see also discussion in EH). As a result we have 
chosen some averaged effective values of mq,mg listed in Table 6 . 

Ea. (j55|l has the form of nonrelativistic Schroedinger equation, however 
one immediately realizes inserting © in (El, that minimization in /ii,p 2 
yields 


min 

Ml 1/^2 


^ ^ /^1 +/^2 

2pi 2/i2 2 2/i J 


-|- mf -|- 


-|- m 


2 

2 , 


i.e. one obtains the Salpeter equation, which accounts for the relativistic 
kinematics. From this simple exercise (valid for zero orbital momentum) 
one can recognize the physical meaning of einbein variables /ri,/i 2 , namely 
Pi ~ (a/p^ + 'rnf) which play the role of constituent masses of the particles 
(for more discussion see jnSlEl)- 

Tn oaTticiilaT Ens.(|54l55() have been used in gni to calculate glueball 
masses at T = 0 with = ar, (which appeared in perfect agreement 

with lattice data EH) where fi for gluons in the lowest mass glueball is cal¬ 
culated to be p = 0.528 GeV. 

Hence one is justihed to apply Eqs. (I54I55|1 both to heavy quark and light 
quark and gluon systems listed in Table 6 . 


6 Results and discussion 


Eqs. (Pi). dsni) have been solved numerically for the parameters shown in 
Table 6 . 

First of all we have concentrated on the charmonium bound states for the 
potential Vi{I) and considered four combinations of parameters Mq, 7 for the 
analytic potential (I48ll5dj) : 


1) Mo = 1 GeV , 7 = 0.69 GeV 

2) Mo = 1 GeV , 7 = 0.2 GeV 

3) Mo = 0.69 GeV , 7 = 0.69 GeV 

4) Mo = 0.69 GeV , 7 = 0.2 GeV. (56) 
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Results of numerical calculations of bound states masses M{cc) for the 
cases l)-4) are presented in Table VII and Fig. 3. As expected, the cases 3) 
and 4) give the deepest bound states which keep intact to largest tempera¬ 
tures, T = 1.6Tc and T = l.STc respectively. For T = Tc charmonium L = 0 
masses he in the interval 3.3 -h 3.4 GeV. 

To follow in more detail the process of bound state dissolution with in¬ 
creasing temperature, we plot in Fig.4 the values of binding energies e = 
M(cc) — 2mc — Vi(oo,T). One can see in Fig.4 that all 4 bound state levels 
are diving in the barrier Vi(oo,T) at some critical value of T/Tc, where the 
radius of bound states is inhnitely increasing, as one can see from the Table 
VII. 

A similar situation occurs for the potential derived from the lattice 

measurements of D^[x) EZII2H1- The masses M[cc) are presented in Table 
VIII for two values of gluon screening mass 7,7 = Mq = 1.044 GeV and 7 = 
0.2 GeV. One can see, that the cc bound states dissolve at lower temperatures, 
1.007Tc and 1.065Tc respectively which reflects the fact, that the amplitude 
of Vi(cxD, T) drops much faster with temperature for the potential Vi{II) than 
for the analytic potential Vi{I). 

Results of numerical calculations for the bottomonium bound state mass 
M{bb) are presented in Fig.5 for the parameter combinations 2 and 4 from 
(jini)- As expected, bound states exist in a wider temperature region up to 
T = 2.2Tc. 

At this point it is possible to compare M{cc) as function of T to experi¬ 
mental masses at T = 0 and to the MEM lattice calculations. 

One can see in Tables VII and VIII that immediately above Tc the values 
of L = 0, M{cc) are about 0.3 GeV higher than the J/'ip mass of 3.1 GeV and 
at higher temperatures M{cc) drop below this level. One can see the same 
type of qualitative behaviour in MEM calculations, (e.g. in Fig. 3 of m and 
Fig.15 of ^2]); however quantitative accuracy of MEM masses is difficult to 
ascertain. 

Of special interest is the L = l,cc bound state. In our calculations for 
the most favourable case 4 this bound state appears at the threshold only for 
the interaction increased 2.5 times (at T = Tc), with the mass around 4 GeV. 
This large difference is explained by the small radius I/Mq ~ 1 GeV“^ of the 
potential Vi(r, T). Thus it is not possible to have P state cc bound states with 
the interaction Vi(/) (and even more so for Vi{II)). At T = 0 the Xco state 
is only 350 MeV above J/t/, which corresponds to larger radius of conhning 
interaction (the radius of Xco is around 0.6 fm). In [T21I221 the mass of Xco for 
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T > Tf. appears below 3 GeV as wide bumps in the MEM analysis. If this can 
be considered as L = 1, cc bound states, the very fact that they appear below 
L = 0, cc states is difficult to explain. The mechanism suggested in jT], which 
takes into account the Thomas spin-orbit interaction, Vs){r) = — ^^^ 2 ^ ■, with 
spacial string tension cth = ^ f D^(x)d^x, cth > cte = 0.18 GeV^ for all 
T > Tc, yields attractive bound states for J = L-|-1 with accumulation point 
just below 2mc = 2.8 GeV. (Note however that this mechanism acts for Xc 2 
rather than Xco)- Let us now turn to other states from the Table VI. 

The interesting light quarkonia L = 0 state with thermal quark masses 
of 0.4 GeV appears to be bound only if one multiplies potential Vi(r, T) by 
the factor 1.33 for the most favorable parameter set 4. The resulting to¬ 
tal mass, as seen from Table 9, is around 1.7 GeV. This is agreement with 
lattice MEM data [SI E] , however one should have in mind, that the ap¬ 
proximate degeneration of S, V, A, P quarkonia states found there, requires 
the use of another relativistic formalism accounting for restoration of chiral 
invariance, and will be reported elsewhere. As it is, the existence of bound 
light quarkonia in 01111 implies that our potentials Vi{I),Vi{II) underes¬ 
timate the actual interaction and should be possibly increased by 30-50%. 
This agrees with derivation of Eq. m, where only the lowest mass Mq term 
is kept, while the next, radial excited gluelump term increases Vi(r, T) by 
roughly 50%. 

Another white state is glueball {gg)i in Fig. 6 and Table 9. One can 
see that its mass is around 2.2 GeV for T = Tc, which is close to the spin- 
averaged T = 0, L = 0 level [101 01]. Lattice glueball MEM data exist 
for T < Tc and are very welcome for T > Tc. Of special interest are 
the colored bound states shown in Fig. 6 and Table 9. One can see that 
the triplet (cc), (eg) and octet glueball survive above Tc- In particular the 
binding energy of {cg)^ state is small, which means that each c quark is 
surrounded by a widespread cloud of gluons, which increases its entropy and 
suppresses the quark Polyakov line. 

All data in Fig. 6 and Table 9 are obtained assuming the potential Vi(/). 
For the potential Vi{II) results are similar, but the bound states disappear 
at smaller T, and in some cases they do not exist for the nominal amplitude 
of interaction. 

Our (cc) and (bb) bound states found from the theoretical and lattice held 
correlators shown in Figs. 3-5 and Tables VII,VIII and IX can be compared 
to the calculations based on phenomenological potentials j^. In 03] 
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the potential was obtained identifying Vi(r,T) with the free energy Fi{r,T), 
measured on the lattice and the resulting binding energies £ for (cc), {bb) 
systems are close to our results shown in Tables VII, VIII and Figs. 3-5. 

Alberico et ah |13j have htted the QQ potential to the lattice free energy 
data and found cc and bb bound states in the interval Tc < T < T^, Td/Tc ~ 2 
for J/'ijj, which roughly agrees with our results, whereas our interval for T is 
much shorter. 

The behaviour of glueball masses is important both from the point of 
view of thermodynamics of the phase transition El Hi, and for the role of 
glueballs in the hadronizing quark-gluon plasma M- 

The hnal comment is on the thermodynamics in the quark-gluon plasma, 
taking into account strong interaction and possible white and colored bound 
states. An extensive discussion of this topic was done in SHl-Eni- In par¬ 
ticular, it was proposed in EH ill, that abundance of colored bound states 
appearing due to strong Coulomb potential (the plateau’s in Vi(r, T) and 
the term aiVi{oo,T) for colored states were not discussed in EH ESI) 
drastically change the plasma thermodynamics. 

However, the presence of self-energy parts for colored states in the form 
of the term aiVi{oo,T) in the total mass (I54|l . can modify contribution of 
colored states to the plasma partition function. E.g., each state appears 
with additional Boltzmann factor exp ^^d for {cg)^ state this 

amounts to the coefficient ~ 0.01 for T c:^Tc = 0.17 GeV. 

Therefore the effective number of d.o.f. in plasma may be much smaller, 
than one could expect. This conclusion is seemingly in agreement with recent 
studies eh "ESI- 
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TABLE CAPTIONS 


Tab. I. Results obtained by fitting simultaneously all the data for the perpen¬ 
dicular electric correlator [see Eq. (P|)] at temperatures T > 
with the functions (EH), where only the non-perturbative coefficient B 
is considered to be temperature dependent. 

Tab. II. Results obtained by a best £t similar to that of Table I, but fixing the 
value of the perturbative coefficient a -|- 6 to the value 0.90/Ac = 0.30 
(marked with an asterisk in the Table) found in [77j . 

Tab. III. Results obtained by fitting simultaneously all the data for the differ¬ 
ence between the perpendicular electric correlator and the parallel 
electric correlator [see Eqs. (ED) and (EH] temperatures T > Tc 
with the functions EH), where only the non-perturbative coefficient B 
is considered to be temperature dependent. 


Tab. IV. Results obtained by a best fit similar to that of Table III, but fixing 
the value of the perturbative coefficient b to the value 0.35/Ac ~ 0.12 
(marked with an asterisk in the Table) found in p7j . 


Tab. V. Parameters of the static potential of binary systems {QQ)d, {QQ)d, 
{Qg)D and {gg)D in different color representations D. 

V^^(r, T) = aV^^{oo, T) + T). 


Tab. VI. Parameters of binary systems, used to calculate bound state energies 
with potentials Vi{I) and Vi{II). Masses are given in GeV. 

Tab. VII. Masses of cc bound states for the potential Vi{I) (in GeV) for four 
different sets of Mo ,7 as in Eq. (EH as function of T/Tc - second 
column. Mean radii of the bound states in GeV“^ ~ third column. 
Binding energies of bound states (in GeV) with respect to the barrier 
height Vi(oo, T), e = M — 2mc — Vi(oo, T) - fourth column. 


Tab. VIII. The same as in Tab. VII, but for the potential Vi(//). 


Tab. IX. The same as in Tab. VII for the systems (ss)i, ( 00 ) 3 , {cg)^, {gg)i, {gg)% 
with parameters given in Tab. VI. 
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Table I 


/9 

T/T, 

B 

M 

(MeV) 

a + b 

XVN 




1042(91) 

0.29(4) 

0.28 

5.90 

1.007 

1.09(34) X 10« 




5.92 

1.030 

0.62(20) X 10^ 




5.95 

1.065 

0.41(14) X 10^ 




6.00 

1.127 

0.18(12) X 10^ 




6.10 

1.261 

0.00(44) X 10^ 





Table II 


P 

T/T, 

B 

(MeV^) 

M 

(MeV) 

a + b 

XVN 




1043(100) 

0.30(>k) 

0.39 

5.90 

1.007 

0.99(33) X 10« 




5.92 

1.030 

0.53(18) X 10^ 




5.95 

1.065 

0.32(11) X 10® 




6.00 

1.127 

0.09(9) X 10® 




6.10 

1.261 

0.00(3) X 10® 
























































Table III 


/9 

T/T, 

B 

M 

(MeV) 

b 

XVN 




887(830) 

0.14(4) 

0.33 

5.90 

1.007 

0.13(49) X 10« 




5.92 

1.030 

0.15(51) X 10^ 




5.95 

1.065 

0.11(44) X 10^ 




6.00 

1.127 

0.06(33) X 10^ 




6.10 

1.261 

0.00(28) X 10« 





Table IV 


P 

T/T, 

B 

(MeV^) 

M 

(MeV) 

b 

XVN 




1154(186) 

0 .12(*) 

0.33 

5.90 

1.007 

0.52(25) X 10« 




5.92 

1.030 

0.54(25) X 10^ 




5.95 

1.065 

0.47(21) X 10® 




6.00 

1.127 

0.37(16) X 10® 




6.10 

1.261 

0.30(12) X 10® 
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Table V 


system 

{QQ)i 

{QQ)8 

{QQ)3 

{QQ)6 

CO 

(Qfi')6 

(^^)l 

00 

CL 

0 

9 

1 

5 

1 

5 

0 

9 



8 

2 

4 

2 

4 


8 

b 

1 

1 

1 

1 

9 

3 

9 

1 



8 

2 

4 

8 

8 

4 



Table VI 


i 

1 

2 

3 

4 

5 

6 

7 

8 

contents 

(&&)i 

(cc)i 

(ss)i 

■ 

icg)3 

(C^)6 


00 

(0 

m\ 

4.8 

1.4 

0.4 

1.4 

1.4 

1.4 

0.5 

0.5 

(0 

m2 

4.8 

1.4 

0.4 

1.4 

0.5 

0.5 

0.5 

0.5 

■ 

0 

0 

0 

1 

2 

1 

2 

5 

4 

0 

9 

8 

h 

1 

1 

1 

1 

2 

9 

8 

3 

8 

9 

4 

1 
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Table VII 



T/Tc 

M 

< r > 


1 ) 

Mo = 1.0 

7 = 



1.0 

3.33 

6.62 



1.1 

3.28 

11.16 


2 ) 

ill 

'o = 1. 

D 7 = 

0.2 


1.0 

3.30 

3.78 

-0.031 


1.1 

3.26 

4.41 



1.2 

3.21 

5.40 

-0.012 


1.3 

3.17 

7.14 

nniiBB 


1.4 

3.12 

11.11 



1.5 

3.08 

23.07 

-0.0002 

3 ) 

Mo 

= 0.6 

3 7 = 

0.69 


1.0 

3.40 

2.77 

-0.107 


1.1 

3.35 

3.06 

-0.077 


1.2 

3.30 

3.45 

-0.053 


1.3 

3.25 

4.02 

-0.033 


1.4 

3.20 

4.91 

-0.018 


1.5 

3.15 

6.67 

-0.008 


1.6 

3.09 

12.13 

-0.002 

4 ) 

Mo = 0.69 7 = 



1.0 

3.36 

2.48 



1.1 

3.31 

2.69 

-0.115 


1.2 

3.27 

2.95 

-0.086 


1.3 

3.22 

3.29 

-0.061 


1.4 

3.18 

3.77 



1.5 

3.13 

4.47 

-0.024 


1.6 

3.08 

5.67 

-0.012 


1.7 

3.03 

8.30 



1.8 

2.98 

18.53 
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Table VIII 


T/Tc M 

< r > 

e 


= 1.0' 

3.39 

3.15 

3.03 

44 7 = 

3.38 

7.08 

20.74 

= 0.2 

-0.041 

-0.006 

-0.0003 
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Table IX 


T/T, 

M 

< R> 

e 

{ss)i 

; Mo 

= 0.69,7 

= 0.20 

1.00 

1.74 

9.76 

-0.0056 

1.05 

1.69 

12.33 

-0.0015 

(CC)3 



= 0.20 

1.00 

3.49 

5.571 

-0.014 

1.05 

3.46 

6.161 

-0.010 

1.10 

3.42 

6.931 

-0.008 

1.15 

3.39 

7.96 

-0.005 

(CC)3 



= 0.69 

1.00 

3.50 


-0.006 

1.05 

3.46 

8.66 

-0.004 

1.10 

3.43 


-0.002 

1.15 

3.39 

12.44 

-0.001 

(C5)3 




■min 


3.55 



2.97 

3.81 


Inra 

2.92 

4.13 


1.15 

2.87 

4.56 


{cgh 

; Mo 

= 0.69,7 

= 0.69 

1.00 

3.06 

4.35 

-0.039 

1.05 

3.00 

4.80 

-0.029 

1.10 

2.95 

5.45 

-0.021 

1.15 

2.89 

6.40 

-0.014 


= 1.0,7 


■min 

2.16 

3.72 



2.11 

4.29 


■nra 


5.12 


1.15 

BB 

6.41 


(55)1 

; Mo 

= 0.69,7 

= 0.20 

1.00 

2.26 

2.10 

-0.335 

1.05 

2.22 

2.21 

-0.289 

1.10 

2.17 

2.34 

-0.246 

1.15 

2.13 

2.49 

-0.207 

{gg )8 

;; Mo 

= 0.69,7 

= 0.20 

1.00 

2.50 

11.27 

-0.003 

1.05 

2.42 

13.48 

-0.0005 
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FIGURE CAPTIONS 


Fig. 1 The behaviour obtained for the static potential Ui(-R, T), given by Eqs. 
(ESI)-®, as a function of the distance R, for different values of the 
temperature T > T^, using for B{T) and M the central values reported 
in Table I and for the parameter b the perturbative estimate ®, with 
the hare lattice coupling constant as = 3/(27r/3) ~ 0.08, as explained in 
section 4. On the right-hand side the scale is multiplied by 7.5, which 
corresponds to ~ 0.6, the characteristic value at large distances, as 
explained in section 5. 

Fig. 2 The three-quark potential Vi{QQ){r,T) in different color states for 
quarks placed in vertices of equilateral triangle with sides r, vs r in 
GeV-h 

Fig. 3 Masses of cc color singlet bound states (in GeV) as functions of T/Tc 
for sets of parameters 1-4 from Eq. dsni). 

Fig. 4 Binding energies of cc,L = 0 bound state, dehned as e = M{cc) — 
2mc — Vi(oo,T) as functions of T/R for the parameter sets 1-4, Eq. 

(jsni)- 

Fig. 5 Bottomonium L = 0 masses (in GeV) vs T/Tc for parameter sets 2,4, 
Eq. dnni), mh = 4.8 GeV. 

Fig. 6 Masses of bound states (in GeV) for the systems labelled in Table VI 
and for the potential Vi{I) vs T/T^ with the parameter set 4 from 
Eq. (®) . 
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